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Abstract 

We present a theoretical calculation on the time- variability of the fine-structure constant to 
fit the result of the recent precise analysis of the measurement of the QSO absorption lines. 
We find the parameters of the scalar-tensor theory to be determined much more accurately 
than fitting the accelerating universe itself, but leading not to easy detections of the effect 
on the equation of state of the dark energy in the earlier epochs. 

To understand the accelerating universe [T] , we have developed a cosmological model in the presence 
of a cosmological constant [5] based on the scalar-tensor theory (STT) in its simplest version. The 
theoretical analysis consists of the two steps. 

In Step I, we derive solutions on the spatially flat Friedmann-Robertson- Walker universe also with 
the assumed spatially uniform scalar field </>, featuring no inflationary expansion in spite of a "large" 
dimensionful constant A which we simply added to STT in its Jordan conformal frame (JCF). As one 
of the most remarkable achievements at this stage, we have successfully implemented the "scenario of 
the decaying cosmological constant," allowing us to understand the well-known number 10~ 120 simply 
as t^ 2 , where to is the present age of the universe, ~ 10 60 in the reduced Planckian unit system with 
c = Ti = Mp(= (87rGo) _1 ^ 2 ) = 1 with Go for today's value of Newton's constant. This conclusion is 
reached in the Einstein conformal frame (ECF) which is identified as the physical conformal frame. 
The scalar field which is canonical here is denoted by a, behaving as dark energy. As we also emphasize, 
the exponential potential of a, a favored choice in the phcnomenological approach of quintessence [3J, 
is an automatic outcome from our theoretically constrained model. 

There are two issues to be discussed. First we argue that the masses of matter fields should 
be time-independent in the physical conformal frame, requiring the traditional Brans-Dicke "model" 
in JCF be replaced by another model, now featuring global scale-invariance except for the A term. 
Secondly, the density of a falls off in the same way as the matter density. This "scaling behavior" 
must be replaced by the "tracking behavior" in order for the universe to be in fact accelerated. This 
situation is realized if the cr-energy stays nearly constant for some duration of time. 

We implemented this in Step II by introducing another scalar field x, which provides a potential 
to which a is temporarily trapped. The field a goes down oscillating toward one of the minima of the 
potential. This oscillation is so small that it hardly shows itself in the process of acceleration, but it 
does in the temporal variation of the fine-structure "constant," through the relation in ECF |2J; 

— = Z—(Aa n 4.6 x 10- 3 A(j, (1) 
a 4tt 



1 



with A always understood for a difference from today's value. Also Z = 5 is for the effective number 
of unit charge of the fundamental quarks and leptons in the loop, while £(= 6 + Alo) = 1.5823 was 
used for our fitting the acceleration. To avoid a conflict with the result of solar-system experiments, 
C 10 _3 [1], we assume a finite force-range of a. We halved the numerical coefficient in (6.194) of 
[2], because we had mistakenly included contributions from the both sides of the photon self-energy 
diagram. Since Aa/a is proportional to the change of a, the former also oscillates, inheriting the 
crucial physical process in the cosmological acceleration. 

The above equation (fT]) is unique in the sense that we derived it based on the standard minimal 
electromagnetic interaction of charged fields. We did not make any ad hoc assumption like the coupling 
in the form of F^F^F^ [5], based on which many phenomenological analyses have been attempted 
on the time- variability of a [3 [7] • We point out that we used the quantum- anomaly technique in the 
relativistic quantum field theory to implement the idea of spontaneous breaking of scale invariance in 
JCF as mentioned before. 

Recently Levshakov et al reported a nonzero time-dependence of the fine-structure constant a 
after a careful analysis of the observed QSO absorption lines, based on the method of Single Ion 
Differential a Measurement; {Aa/a) x 10 6 = -0.12 ± 1.79 and 5.66 ± 2.67 at z = 1.15 and 1.84, 
respectively [51 O [TU] , corrected also for the improved estimate of the Fe II sensitivity coefficients 
[i~T] . We note that the quoted errors include both statistical and systematic uncertainties. Taking a 
rather tentative nature of their analysis into account, we still believe it worth attempting to apply 
a theoretical estimate (fTJ) to their result, going beyond our earlier phenomenological analyses [T2] on 
the then available observational results [13] . We also include the constraint obtained from the Oklo 
natural reactors [14j as well as the laboratory measurement of a/ a [15] . In this way we expect to 
constrain the parameters and initial values much better than fitting the accelerating universe itself. It 
seems unlikely, on the other hand, that this type of solution results in the past-time deviation of the 
equation of state w = p/p for the dark energy from —1 at any practical level of detection, indicating 
probably an important insight into what behavior of w is going to look like. 

As we admit, however, we have never swept the entire space of 5 parameters and 4 initial values, 
having shown the fit in Fig. 5.8 of [5] only to illustrate an explicit solution that reasonably fits the 
data for the accelerating universe. This solution, to be referred to as the "reference" solution, turns 
out, however, to predict too large values on Aa/a to be compared with the recent observation [8"]-|10|. 
As a remedy we change the value of the parameter 7, a coefficient in front of the sine-Gordon term in 
the potential in (5.58) of [2J, from 0.8 to 0.81. Corresponding changes of other parameters and initial 
values are to be in order. In view of the nonlinear nature of our cosmological equations, however, we 
follow a cautious approach by starting with the reference solution, with the parameters (A — 1, £ = 
1.5823, to = 4.75,7 = 0.8, k = 10.0) and the initial values (d = 6.7544, a\ = 0, xi = 0.21, xi = 0) 
imposed at t\ = 10 10 (sometime after the primordial inflation), also in the reduced Planckian unit 
system. We search only the neighborhood of the reference solution except for 7. 
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The solutions we reached in this heuristic manner are shown in Fig. 1, with two fixed parameters 
p and Sx and the initial value a\ varied in the range including the reference value 6.7544. Let us begin 
with short explanations on the proposed constants p and Sx. 

We are going to identify the epoch about 1.9 x 10 9 years ago for the Oklo phenomenon with a 
zero of the oscillating Act /a. This is affected obviously by a factor p chosen to give the age of the 
universe to = p x 1.37 x 10 10 y. We also notice that the calculated curves for Aa/a will be shifted 
horizontally in Fig. 1, if we change the initial time t\. More technically, our cosmological equations 
were integrated with respect to x = log 10 t starting at the initial value x\ chosen conveniently to be 
10. By changing x\ by a small amount Sx, we are allowed to re-use the result of the reference solution 
but with 7 = 0.81, resulting in the approximate shifting the curves in proportion to 5x. From an 
overall fit to the Oklo at around z sa 0.16 together with two QSO data, we favor the choice p = 1.02 
and Sx — 0.05 as reasonable fine-tunings. The former choice for to » 1.40 x 10 10 y will be tolerated. 



s 
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Figure 1: (Aa/a) x 10 6 is plotted against the redshift z for different initial values o\ as shown in 
each entry, to be compared with the QSO data [5]- [TO]. Also y' = (d(Aa/a)/dz) z= o, the slope of the 
calculated curve at z — 0, is to be compared with the laboratory measurement of a/ a |15j . The Oklo 
constraint [14] is marked by a blob at z rs 0.16. 



Note also that the Oklo phenomenon gives an upper-bound of |Aa/a| < 10~ 7 , or even smaller 
<, 10 -8 [14] . represented by a blob around z w 0.16, with a nearly invisible error bar in this plot. 

Given estimates of p and Sx as above, we now vary o\ , which might be considered as a representative 
of other initial values or parameters. In Fig. 1, we find quite a variety of the behavior of y(z) = Aa/a 
for different values of o\ in a rather short interval of 0.00015 including the reference value 6.7544. 

We also note that the slope of the curve y' at the origin is related to the result of the laboratory 
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measurement of (d/a)o at the present epoch by 



- =-Hoy , a , 



(2) 



as given most recently by (—0.26 ± 0.39) x 10 15 y 1 [T5], which can be translated into y' = (3.5 ± 
5.2) x 10~ 6 for h = 0.73, also to be denoted by y[ ht ± tribt- 

As we find, the computed curves agree reasonably well with the two QSO data at z — 1.15 and 
z = 1.84 for <7i > 6.7544, which happens to be the reference value in the solution (c), but with 
unacceptably large y' Q with v = (y' — Z/i'bt) / °"ibt £4.5. The deviation v is larger for larger a%, while 
it is smaller for smaller a\, coming down to v = 0.3 in (f), for example, but now with an obviously 
poor agreement with the QSO data. In spite of a desperately necessary effort to search for wider class 
of solutions, also given a still small number of precise measurements particularly on the QSO result, 
the very presence of "compromises," exemplified by the solutions (d) or (e), with v = 2.8 and 2.0, 
respectively, seems encouraging since, as we show shortly, all the solutions discussed here share nearly 
the same cosmological consequences. 

More explicitly, we notice that Figs. 2 and 3, corresponding to Figs. 5.8 and 5.9 in [3J, respectively, 
on the cosmological behaviors computed from the typical fit (e) in Fig. 1, remain nearly the same for 
all other solutions (a)-(f); no discernible differences are recognized for their results. This is a general 
tendency, due to the fact that considerable difference in Aa/a comes from the same in a, which is 
only a "small" ripple if viewed from the overall cosmological behaviors, as shown in these plots. 
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Figure 2: Cosmological behaviors in ECF for the 
solution (e) of Figure 1, with the reference coun- 
terpart in Fig. 5.8 of [2J, p and p s for the densities 
of the matter and the a\ system (dark energy), 
respectively. Present epoch, log 10 to — 60.20, is 
shown by a vertical dotted line. 
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Figure 3: Magnified view of Figure 2 near the 
present epoch, corresponding to Fig. 5.9 of [2]. 
Today's values of h and Q\ are also shown, mea- 
sured according to the right-hand vertical scale, 
found within the observed ranges [TB] . 



In fact, corresponding to Fig. 5.10 of [2], we have Fig. 4 with the vertical scale magnified by as 
much as 10,000 times compared with that of the upper panel of Fig. 2. Measuring Aa/a is something 
like watching a portion of a around the present epoch, featuring an apparently smooth and flat section 
followed by a sudden but tiny step-like rise, in the upper panel of Fig. 2, through a microscope with 
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much better resolution. In other words, Aa/a, if confirmed, is expected to probe the parameters 
much more accurately than the cosmological acceleration can do. Just for the sake of comparison, we 
show the corresponding reference solution in Fig. 5, basically the same as Fig. 5.10 of [2]. 
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Figure 4: Magnified view of the upper panel of 
Fig. 2, with the vertical scale magnified by about 
10,000 times, a is being trapped by the potential 
provided by x? which starts falling toward zero, 
thus releasing a eventually. 
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Figure 5: The reference solution counterpart of 
Fig. 4. The behaviors may look similar, but the 
rate of vertical magnification is roughly 30 times 
less than in Fig. 4. 



The difference between Figs. 4 and 5 can be understood by different values of a when it has 
"cruised" before it is captured by the potential due to x- We simply say that a in Fig. 4 "happened" 
to be located much closer to a potential minimum than in Fig. 5. This is the reason why the ripple of 
a is sufficiently small to be in agreement with the small observed Aa/a. In other words the reference 
solution would have predicted Aa/a an order larger than the current observation [8]-[10j. 

Figures 4 and 5 also suggest that simple behaviors of a in the immediate past might have left 
some of observable effects. In principle, the oscillation of a entails a more or less oscillatory behavior 
of w = p s /p s = (K - V)/(K + V), with K = cr 2 /2 + x 2 /2, while V is given by (5.58) of [2J for 
the interaction energy of the system of a and x which together comprise dark energy, denoted by 
the subscript s, in Step II in our theoretical model. Unfortunately, the relationship between the two 
quantities appears to be complicated as suggested by comparing Figs. 1 and 6, partly because x is 
also involved. Furthermore the magnitude of time-variability of w is much smaller than in smooth 
behaviors which are expected on phenomenological bases [7J [UJ [TB] ■ This is a general feature shared 
by other solutions in Fig. 1 as well. We learn that apparently structureless behavior of w close to —1 
in a range of small z may not necessarily imply a purely constant A. It can still be consistent with a 
nonzero Aa/a of the order of 10~ 6 . 

This apparently flat behavior stuck closely to —1 in a limited resolution is shown to continue down 
to t ~ 10 54 at which p s begins to rise in the past direction, as we illustrate in Fig. 7. This epoch 
nearly coincides with the "equal time," or the time toward the end of CMB epoch. Nearly at the 
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same time, the effective power exponent a = tH, allowing a locally approximate dependence a ~ t a , 
begins, in the future direction, to change from 1/2 to 2/3. In fact we adjusted the parameters and 
initial values in such a way that the condition p s <C p is met around t ~ 10 , thus maintaining the 
success of the primordial nucleosynthesis undisturbed by a cosmological constant, as also seen in the 
lower panel of Fig. 2. It may not be easy to change this feature expected from an overall point of 
view. We still find it likely that the smallness of the ripple of a is responsible for near undetectability 
of the w variability, as we inspect the reference counterpart in Fig. 8. 

Redshift z 
6 4 3 2 1 

-0.99 
-0.993 
-0.994 
-0.996 
-0.998 
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log t 

Figure 6: The equation of state w — p s /p s for the solution (e) of Fig. 1 stays so close to the lower 
bound —1 that no detection appears feasible for reasonable range of z. Basically the same result 
follows for any of other solutions displayed in Fig. 1. 

Unlike in Fig. 7, we do find a number of spike-like behaviors of w. The biggest one centered at 
log 10 t ~ 58.95, also including smaller ones toward log 10 i ~ 60, is due to the slight change of p s , as 
also shown in Fig. 8. The same type of behavior of p s does not seem an uncommon occurrence among 
many examples of solution, though rarely occurring for such low values of z < 2, or log 10 t > 59.6. It 
might be rather accidental to find no such spikes, though allowing to be much closer to —1, in the 
solutions shown in Fig. 1. The behaviors of the sharp rise (in the past direction) to w = +1 around 
log 10 t ~ 54 are nearly the same as in Fig. 7. 

Perhaps more important is to note that major behaviors, like the above examples, but unlike those 
discussed in pj I17j from phenomenological points of view, may show the features having something 
in common with a delta function or a step function, to which a power series expansion is known to be 
ineffective. This might be a unique conclusion we draw from our theoretical model designed primarily 
to understand the accelerating universe. We should be prepared, however, at least for the possible 
occurrence of such "singular" terms. 
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Figure 7: Inter-relations among the matter den- 
sity p, the a\ system density p s (dark energy), 
measured according to the left-hand scale, while 
the equation of state for the scalar fields w, and 
the effective power exponent a of the scale factor, 
allowing a ~ t a , measured against the right-hand 
scale, for the solution (e) in Fig. 1. 



Figure 8: The same as in Fig. 7, but for the 
reference solution. Note the occurrence of con- 
siderable time-dependence of w, unlike in Fig. 7. 



As we emphasized before, we are not fully satisfied by the solutions obtained in our limited and 
heuristic approach. Also the future re-analysis of the QSO measurements may reveal aspects different 
from those in [S]. Observational searches on the equation of state available at present, on the other 
hand, appear to be consistent with nearly anything including purely w = — 1. Taking these circum- 
stances into account, we should make further efforts to search for other types of solution beyond the 
neighborhood of the reference solution. 



7=.81, (7, =6. 75431 
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Figure 9: An example of the "desired" behaviors, instead of those in Fig. 1, obtained artificially by 
shifting the time coordinate x in Aa/a from the value responsible for the acceleration by Ax — 0.534, 
agrees reasonably well with all of the observations, including v = 0.4. 



It might be still interesting to suggest that accepting curves as shown in Fig. 1 may not be the 
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only way to fit the available data. In fact Fig. 9 shows an example of the "desired" behavior, which 
we constructed artificially by a mathematical trick to use a(x — Ax) in ([T]) where Ax = 0.534, by 
exploiting a low-frequency portion of the oscillating a. With o\ = 6.75431 we find that the curve 
quite different from those in Fig. 1 achieves reasonable agreements with all the observations including 
the laboratory experiment, though with the calculated w even closer to —1 than in Fig. 6. We expect 
to see how this behavior emerges in fact, by either changing the parameters and initial values in wider 
ranges, or modifying the potential slightly in (5.58) in [2]. 

Searching for solutions beyond the reference solution might be suggested by yet another indication 
of the dark energy density to be present for z >, 2 for the required growth of structure [15] , 
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